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The phase diagram of a stack of tensionless membranes with nonlinear curvature 
energy and vertical harmonic interaction is calculated exactly in a large number of 
"q" . dimensions of configuration space. At low temperatures, the system forms a lamellar 

o 

phase with spontaneously broken translational symmetry in the vertical direction. At 

J>y a critical temperature, the stack disorders vertically in a melting-like transition. The 

critical temperature is determined as a function of the interlayer separation /. 
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I. INTRODUCTION 



Recently, a model for a finite stack of tensionless membranes [[]]] was studied with respect to 



the effects of higher-order terms of the curvature energy [Eh. The approach was perturbative, using 

s 

the renormalization group to sum infinitely many terms. It was shown that thermal fluctuations 



induce the melting of the stack into a vertically disordered phase. By its nature, the perturbative 
expansion was able to give a satisfactory description only for the ordered phase. 



^ ■ For a description of the disordered phase and a better understanding of the entire transition, 



we analyze in this paper the behavior of a stack of tensionless membranes exactly for very large 
dimension d of the embedding space. Since the model is exactly solvable in this limit, we can 
calculate all its relevant properties explicitly, in particular its complete phase diagram as a function 
of the interlayer separation /. 
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II. THE MODEL 



As in Ref. [fZp, we consider a model in which a multilayer system is made up of (N + 1) fluid 
membranes, parallel to the xy-plane of a Cartesian coordinate system, separated a distance /. If 
the vertical displacement of the mth membrane with respect to this reference plane is described by 
a function u m (x) = w(xj_, ml), where xj_ = (x, y), the energy of the stack reads: 

Bn 



f ±J 



(1) 



Here, H m = d a N m ^ is the mean curvature, where N m oc (—diU m , —d 2 u m , 1) is the unit normal 
to the mth membrane, and 

9m,ij &ij diU m OjU m (2) 

the induced metric, with i,j = 1, 2, d\ = d/dx, d 2 = d/dy and g m = det[g m> ij). The parameter 
k is the bending rigidity of a single membrane, and B the compressibility of the stack. In Eq. 
([]]), as in the following, the subscript denotes bare quantities, whereas renormalized parameters 
will carry no subscript. 

For slow spatial variations, the discrete variable ml may be replaced with a continuous one, 
and u(x_l, ml) — > u(x), where x = (xj_, z). In this limit, the energy (jl]) reduces to 

E = dz J d 2 x ± ^ [\K,H 2 + \B,{d z uf] . (3) 

Here we have introduced the bulk version of the bending rigidity K = kq/ I, and defined L\\ = Nl. 
The gradient energy (d z u) 2 has the unphysical feature that it gives ^-dependent reparametrizations 
a kinetic energy. It should therefore be replaced by the normal gradient energy (N ■ Vu) 2 . We have 
seen in Ref. [|2|] that this changes the critical exponent with which the renormalized B vanishes 
as the bending rigidity becomes critical. However, in the limit of large d, which we are going to 
investigate, the difference between the two gradient energies will be negligible. 
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III. LARGE-ci ANALYSIS 



For arbitrary d, the vertical displacement of the mth membrane in the stack becomes a (d — 2)- 
vector field u m (xj_). 

A. Partition function and the energy 

It is useful to consider as an independent field [Q], and impose relation @ with help of a 
Lagrange multiplier Ay. We write the partition function as a functional integral over all possible 
configurations u m (xjj of the individual membranes in the stack, as well as over all possible 
metrics g m ,ij- After taking again the continuum limit, the partition function reads: 

Z = J VgV\Vue- Eo/kBT , (4) 

with 

E = J dzd 2 x x ^g{a + \B {d z u) 2 + ±K Q (d 2 k u) 2 + \K^{5 i3 + c^u - 9ij ) - |t A^} , 

(5) 

where u is a (d — 2) -dimensional vector-function of xj_, z. Note that the functional integral over 
A in @ has to be performed along the imaginary axis to result in a 5-function. We have also 
introduced a term proportional to \ 2 V This term is necessary to renormalize the theory, and its 
coefficient r corresponds to the large-d in-plane compressibility of the membranes. Since we 
take the membranes to be incompressible, we shall set the renormalized r equal to zero at the end 
of our calculations. We have also included a surface tension cr , again to absorb infinities and to 
be set equal to zero after renormalization. 

The functional integral over u in Eq. (Q) is Gaussian and can be carried out to yield an effective 
energy 

E eS = E + E 1 , (6) 

with 
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E = J dzd 2 x±^/g [a + \K Q X ij {8 i5 - g^) - \t q \ 2 



(7) 



and 



Ex 



d-2 



fc B TTrln B u 2 + K (q{ - ft A%- 



(8) 



where the functional trace Tr is here an integral over space as well as the integral over wavevectors 
q_i_ and u, after replacing d 2 — > — lu 2 and g'^didi — > — q 2 . Note that the discrete nature of the stack 
restricts the integral over the wavevectors u to the first Brillouin zone \u\ < it /I. 

In the large-d limit, the partition function (Q) is dominated by the saddle point of the effective 
energy @ with respect to the metric and the Lagrange multiplier A u . For very large mem- 
branes, the saddle point can be assumed to be symmetric and homogeneous 



Qo 



(9) 



with constant ^ an d ^o- At the saddle point the effective energy @ becomes the free energy of 
the system. 

In the following we shall investigate both the case of an infinite and a finite stack of membranes. 
As we will see, the large-d approximation allows for the vertical melting even in an infinite stack, 
which is not found perturbatively [[]]]. 

6. Infinite stack 



Let us first analyze the case of an infinite stack. To simplify our calculations, we assume the 
number iV + 1 of membranes in the stack to be very large, making the distance / between them 
very small. In this regime, we may extend the limits ±n/l of the integral over uj to infinity. The 
explicit /-dependence will be introduced later into our calculations. 

After evaluating the functional trace in Eq. (Rf), we obtain 



jp dk B T r 2 K JA 4 Ao.o 
Ex = — — / dzd x± qqJ—< — + —A + 



\ 2 
A 



2 In 



4AJ 

An 



(10) 



2 J V B Q {8tt 8tt 64tt 

where ultraviolet divergences are regularized by introducing a sharp transverse wavevector cutoff 
A and d — 2 has been replaced by d for large d. 



We may now absorb the first term in < JT0| ) by renormalizing a , so that 

dk B T [Kq . 4 



a = a + 



'—A 

167T V B 



(11) 



is the physical surface tension, which is set equal to zero. The second, quadratically divergent 
term in < JlQ| ) is used to define the critical temperature as 

1 <fe A* 



T c 16tt 

The next divergent term, proportional to Aq, is regularized by introducing a renormalization scale 
fi and modifying the the in-plane compressibility to 



(13) 



32vr f5 \ 

The physical in-plane compressibility r is now set equal to zero, as explained in the previous 
section. 

The effective energy thus becomes 



E eS = Jdz d 2 x ± KXg | - lj + ^ + 



A 



(14) 



with the constants a = g?/cb/647T\/B, A = / u 2 e -1 / 2 . 

From the second derivative matrix of E^ff with respect to £> and A we find that the stability of 
the saddle point is guaranteed only for A < A. 

Extremizing the above expression with respect to g, we find two solutions for A, namely A = 



and A = Aqo, with 



k n i 

f ~ T r 



(15) 



For T < T c , this equation has no solution for X^. In this case, the only possible solution is 
A = 0, which corresponds to the ordered phase as we shall verify later. For T > T c , the saddle 
point lies at A = A OT , which is now well-defined. This is the vertically disordered phase. 

The free energy density at the extremum is given by 



f = K\ c 



(16) 



and its behavior is similar to the one found perturbatively in Ref. [TJ]. (see Fig. [[]). 



/ 




FIG. 1. Free energy density of an infinite stack of membranes. The dashed curve indicates the unphysical 



branch of the solution of Eq. ( |15| ) 



Extremizing the effective energy (|T4|) with respect to A, we find g as a function of temperature. 
For T < T c it is given by 



= 1 - 1- 

T 

± c 



(17) 



This as T approaches T c from below, indicating the vertical melting at T c . In the disordered phase, 
q is found to be 



Q- 



-' = l-i 



aXooT 



(18) 



VK 

As T approaches T c from above, Aoo tends to zero, and g goes again to infinity. 

The positivity of g and the stability of the saddle point imply that there is a maximum temper- 
ature, given by 

— = --— (19) 

± max ± c V 

below which our assumption that the membranes in the stack are in-plane incompressible does not 
lead to a stable system. 



6 



C. Finite stack of many membranes 



Let us now analyze the case of a finite stack of size L\\ . Now the functional trace in i 
a sum over the discrete wavevectors u n , given by 

u n = —n, n = 0, ±1, ±2, • • • . 

For small A , a series expansion leads to 

Ei = 

with 



involves 



(20) 



2 dk B T 
dza a — - — 



£o ei, 



(21) 



A 4 jK^ 7i 
61 ~ 8^V^o ~ 12 V K L 2 



B 1 Aq ^o_p2 _|_ 



Ao ln /^ 

47tLh I Lii 




An #0 



+ 64ttV 5 



3 - 27 + 2 In 



An -t^ 



00 / j^m+l^m 



8ttA 2 Ly V i^o, 
5o 



<(m-l). 



(22) 



m=3 m^<» 7T» " \^ / r(f) 
As in the case of the infinite stack, we absorb the logarithmic divergence by renormalizing the in 
plane compressibility via Eq. (|T3"|), setting r equal to zero for incompressible membranes. The 
surface tension receives now an Lu -dependent renormalization 



O = (To + 



dk B T / Kq 



71 Ba 1 



16n V B 24 V KqL° 



2 > 



(23) 



and a is again set equal to zero to describe a stack of tensionless membranes. 

Extremization of the renormalized combined effective action (0) and ( PH ) with respect to g 
leads again to two possible solutions for the saddle point, namely A = or A = Xu , with 



A, 



In 



+327T 

where 



3/2 



, A 

00 

E 

m=3 



1 

~ 2 
-l) m+1 A m_1 



,L] 1 IB 



m 2 2m n m 



-L] 



m-2 



"i-2 „ 

— r 



m— 1 ' 



r(f; 



-C(m 



v 7 ^ (I 1 
a \T T L|I 



(24) 
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(25) 



is the inverse critical temperature for a stack of size Ly . 

For T < T L . , Eq. ( |24| ) has no solution. In this case, the stack is in the ordered phase, the only 
available solution for the saddle point being A = 0. For T > T L .., there exists a nonzero solution 
\u , where the system is in the vertically disordered phase. 

Let us now examine the saddle point solutions for g. In the vertically disordered phase where 
A = Al,i is nonzero, we may expand the effective energy in a small-Ly series. Extremization with 
respect to X L leads to 

-i T , aA L|| T 



2K V ^ 3 2 2m 7T m V 771/ " \BJ T( 2 , 

The positivity of g and the stability of the saddle point again define a maximal temperature, given 
by 



i i dkB T j^vm (27) 



TmL ^max 167rKL|| V^ B T c A 

above which our assumption that the membranes in the stack are in-plane incompressible cannot 
be maintained. 

In the ordered phase, the situation is more delicate. For A = 0, g can be calculated exactly, and 
we obtain 

smh [ 



2L^ 



(28) 



with an infrared regulator equal to the inverse lateral size of the membranes in the stack. If the 
size Lii of the stack is large, ( p8| ) may be approximated by 

Q- 1 « 1 - ^. (29) 

For smaller stacks, however, the positivity of £ is not guaranteed. For a fixed, but small stack size 
Lii, and for fixed lateral size L± of the membranes in the stack, there is a characteristic temperature 
defined by 

8 



T* = ^S- , (30) 

dk B \n(lGTiVBKLi/dk B T c ) 



above which g changes sign, and ( |28| ) is no longer applicable. Interestingly, for all L± and for 
all finite sizes L\\ of the stack, the critical temperature T L is lower than T*, so that the vertical 
melting still occurs. The behavior of g is depicted in Fig. ^. 




T r „, T* T r 



FIG. 2. Behavior of g as a function of T. The solid lines indicate the solutions of the saddle point for 
q^ 1 for a finite stack. Above T^,,, is given by (26), and below T^,, by (28). The dashed lines indicate 



the behavior of g for an infinite stack. 

Note that Eq. ( PD] ) reflects the existence of a characteristic horizontal length scale. At fixed 
temperature T L <T <T*, and for membranes of lateral size L± smaller than 

the height fluctuations of the individual membranes are not strong enough to destroy the ordered 
phase. The characteristic length L p corresponds to the the de Gennes-Taupin persistence length £ p 
[0] of the individual membranes, below which crumpled membranes appear flat. 



D. Finite number of membranes 

Until now we have performed our calculations in the somewhat unphysical continuum approx- 
imation, by letting the interlayer separation / be very small making the number of membranes in 
the stack very large. Let us now investigate the properties of the stack for a fixed number + 1 
of membranes at a finite interlayer distance /. 
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For this purpose, we replace the continuum derivative d\ in the z-direction with the discrete 
gradient operator V 2 , whose eigenvalues are given by 



V 2 s 



2(1 — cosu; n Z) 
P 



where g is some test function. The discrete wavevectors uj n are now given by 

nit 

Un = WV n= 1 ' 2 '"'' 

For small interlayer separation I , the free energy is given by ( |2"T| ) with 



A 4 K 1 



A|_ Ko 
647rV B 



Bp Aq Kq 2 Aq 



lniV + 2iVm UW ^A 2 



2 In 



'4A 2 

An 



i v (-i) m+1 A^ /m _ 2 /^o^ r(^' 



(32) 



(33) 



r(f) 



Gv(™ - 1), 
(34) 



where we have defined the modified Zeta-function 



N 



n=l 



(35) 



I 1 - «* (t ; 

We proceed by renormalizing the in-plane compressibility via Eq. (|T3|), setting r equal to zero 
for incompressible membranes as before. The surface tension receives an Z-dependent renormal- 
ization 



dk B T [Kq 4 d/c B T [Bq 



(36) 



and a is set equal to zero to describe a stack of tensionless membranes, as before. But now the 
bulk bending rigidity is also modified to 



Note that this renormalization agrees with the known result for a single membrane 
The saddle point for A is now given by A = or A = A/, with 

AA 1 



(37) 



A, 



In 



a 



X) 2 

K ( 1 1 
T ~ % 



oo / i\m+l\m-l 



m=2 



- 3m+l 

ml 2 



~BJ T 



r(2 ti ; 



CAr(m 



(38) 
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where 



1 



1 dk-Q 
¥ r + 8itNk 



— In TV + 2N 




(39) 



is the inverse /-dependent critical temperature, and n is the bending rigidity of a single membrane 
in the stack. The two solutions for A again imply the existence of two different phases, with a 
phase transition at the critical temperature T\, which, as in the perturbative case, depends only 
weakly on the number of membranes in the stack. The corresponding solutions for g, obtained by 
extremizing the effective energy with respect to A, are given by 



for T <Ti, that is, in the ordered phase, and 

aXiT 



Ti 



(40) 



T 

Ti 



1 



A' 

oo 



dk B T ~ (-l) m+1 Ar _1 



2 \ , m -2 



1--W 

mj 



b) T(fT Cjv( 



(41) 



in the vertically molten phase, where T > 7]. 

The stability of the saddle point requires a minimum interlayer separation 



It, 



> / B , T j_ / Atxk 
■ ' — N 2N exp 1 



(42) 



K ^ \dksTc 

below which the stack becomes unstable. l m ™ is inversely proportional to de Genne's penetration 



depth y K/B [Q], which is of the order of the interlayer separation in smectic liquid crystals. 
The phase diagram of the stack is depicted in Fig. |3] 




Disordered phase 



Stack 



FIG. 3. Qualitative phase diagram in the I x T plane. The critical line is plotted for I > l m in, for a fixed 
number N + 1 of membranes in the stack. As I increases, the critical temperature T\ goes asymptotically to 



zero. 
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E. Properties of the phases 



Let us now characterize both phases in more detail. For temperatures higher than TJ, the 
solution of the saddle point is A = A;. This corresponds to the disordered phase, where the stack 

— 1/2 

melts. The normals to the membranes are uncorrelated beyond a length scale X t , as can be 
derived from the expression for the orientational correlation function, which in the limit N — > 
oo, I — > 0, with constant Nl = L\\, reads: 

(0 f u(xj., z)d jU (x'±, z)) ~ %e-^ x ^- x 'xl. (43) 

The length scale A ; may thus be identified with the persistence length £ p [p~0t] - 

In the low-temperature phase, the solution of the saddle point is A = 0. This corresponds to 
the ordered lamellar phase as can be seen by examining the orientational correlation function in 
the planes of the membranes. We find in the limit iV — > oo, I — > 0, with constant Nl = Lf 

(d iU (x±, z)d jU (x'±, z)) ~ i 5 ^ Tr (44) 

l X -L x ±l 

This slow, algebraic fall-off of the correlation function implies that, at large distances, the normal 
vectors to the membranes are still parallel, so that the surfaces remain flat on the average. The 
effect of thermal fluctuations is suppressed, and they do not disorder the stack. 

Let us now calculate the entropy loss in the ordered phase. By a simple scaling argument [jlp, 
this quantity should be inversely proportional to the quadratic interlayer spacing, 

-TAS = ^ } (45) 

where a is a temperature dependent proportionality constant. As we shall see, as the interlayer 
distance increases, logarithmic corrections must be added to (|4"5|). 

The entropy loss in the ordered phase can be computed by calculating the difference between 
the free energy density of a single, isolated membrane, an the free energy density of the stack. For 
small values of A, it is given by 

1 l~B / 7T \ 1 A A , AI [K\ A * . nn , d „ 

" TAS = 5V if ( + COt *N) W ~ M i 1 - TV 1? j + 2^S" 1Sm IN' <46) 
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Strictly speaking, the ordered phase corresponds to A = 0. In that case, ( pq ) agrees with (ftq), 
and we see no correction to the entropy loss. However, if the size of the membranes in the stack 
is smaller than the characteristic length L p , an ordered phase still exists for small values of A [see 
discussion after Eq. (Plf)], in which case the corrections to the first term in ( |4"6] ) will appear. 



IV. CONCLUSIONS 

In the limit of large embedding dimension d, we have shown that a stack of tensionless and 
incompressible membranes melts vertically upon approaching a critical temperature, where the 
lamellar phase goes over into a disordered phase. In contrast to the low-temperature ordered phase, 
where the decay of orientational correlations is power-like, the high-temperature disordered phase 
is characterized by an exponential decay of orientational correlations, with different length scales 
in the transversal and longitudinal directions. 
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